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CN I The study of high energy nuclear collisions has entered a new stage with RHIC; it 

Qh| therefore seems a good time to ask what we have learned from the experimental results 
• obtained up to now. I recall what we had expected to find when the SPS and AGS 
programs were started, summarize what actually was found, and then try to assess what 
we have learned from the results. 
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! 1. The Beginning 

I , It all began with the idea of an intrinsic limit to hadron thermodynamics. During the 

0\ \ past fifty years, different conceptual approaches had led to an ultimate temperature of 
^ ' strongly interacting matter. Pomeranchuk [ |l]] first obtained it from the finite spatial 

O ■ extension of hadrons: a hadron can only have an independent existence if it has an 
^ . independent volume. Then Hagedorn [|^ arrived at a limiting temperature by postulating 
Q-i! a self-similar hadronic resonance composition: a resonance consists of resonances which 
Qh! consist of resonances, and so on. The resulting excitation spectrum was later also derived 
in the dual resonance model [ |^ . With the advent of the quark infrastructure of hadrons 
and of quantum chromodynamics, it became clear that the ultimate temperature found 
in all these considerations was really a transition point to a new state of matter, to a 
plasma of deconfined quarks and gluons [ ^ . 
5^ \ Statistical QCD, in particular in the finite temperature lattice formulation, has sub- 

sequently confirmed this hypothesis: at a now rather well determined temperature (for 
vanishing baryon density, ~ 150 — 180 MeV), strongly interacting matter undergoes 
a transition from a medium of color singlet hadronic constituents to one of deconfined 
colored quarks and gluons [ |^. The energy density at the transition point was found to 
be e(Tc) ~ 0.7 — 1.0 GeV/fm^. Moreover, the transition turns a hadronic state of spon- 
taneously broken chiral symmetry into a quark-gluon plasma in which this symmetry is 
restored: at Tc, the effective constituent quark mass of some 0.3 GeV vanishes, and we 
recover the bare quark mass of the QCD Lagrangian. 

The obvious desire to test this fascinating phase structure of strongly interacting matter 
first led to the fixed target experiments at the AGS in Brookhaven (with y/s ~ 5 GeV) 
and at the CERN-SPS (with ~ 20 GeV). In 1986/87, light ion beams on heavy ion 
targets started the program, and in 1994/95, heavy ion beams followed. Today, much of 
this program is concluded. So, what have we learned during the past fifteen years? In 
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this opening talk, I will address that question by asking: 

• What did we expect to find? 

• What did we find? 

• What does that tell us? 

In my report, I will first recall briefiy the expectations concerning signatures at the be- 
ginning of the experimental heavy ion program at the AGS and SPS in 1986 and then 
summarize what had really been found when it was (in first generation experiments) com- 
pleted in 2000. Following this, I will try to indicate what conclusions can be drawn from 
these results, for the conditions reached, from the hard probes of the early stages and 
from the observed hadronisation pattern at freeze-out. 



2. Hopes 

The evolution of a high energy nucleus-nucleus collision was pictured in the form shown 
in Fig. |l]. After a rather short equilibration time tq ~ 1 fm, the presence of a thermalized 
medium was assumed, and for sufficiently high initial energy densities, this medium would 
be in the quark-gluon plasma phase. 




The initial energy density of the produced medium at the time of thermalization was 
(and still is) generally determined by the Bjorken estimate [ ^ 

- ( Wh 

^ V '^y )y=o ^RW 

where {dNh/dy)y=o specifies the number of hadronic secondaries emitted per unit rapidity 
at mid-rapidity and Wh their average energy. The effective initial volume is determined in 
the transverse plane by the nuclear radius Ra, and longitudinally by the formation time 
To of the thermal medium. 

The temperature of the produced medium was assumed to be observable through the 
transverse mass spectrum of thermal dileptons and the momentum spectrum of thermal 
photons [ ^ |[. The observation of thermal dilepton/photon spectra would also indicate 
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that the medium was indeed in thermal equilibrium. The functional form of the spectra is 
the same for radiation from hadronic matter and from a QGP; but the observed rates and 
temperatures were expected to differ in the two cases. It was clear from the beginning 
that these signals would be emitted during the entire thermal evolution of the system, 
making a separation of different phases of origin not very straight-forward. 

The determination of the nature of the hot initial phase required a signature sensitive to 
deconfinement. It was argued that in a deconfined medium the ^/ip would melt through 
color screening [ ^ and that, therefore, QGP production should lead to a suppression of 
J/'0 production in nuclear collisions, compared to the rates extrapolated from pp data. 
Similarly, the QGP was expected to result in a higher energy loss for a fast passing 
color charge than a hadronic medium, so that jet quenching [ |lOl| should also signal 
deconfinement. 

The behavior of sufficiently short-lived resonances, in particular the dilepton decay of 
the p, was considered as a viable tool to study the hadronic medium in its interacting 



stage and thus provide information on the approach to chiral symmetry restoration [ |TT 

The expansion of the hot medium was thought to be measurable through broadening 
and azimuthal anisotropics of hadronic transverse momentum spectra [ The size and 



age of the source at freeze-out was assumed to be obtainable through Hanbury-Brown- 



Twiss (HBT) interferometry based on two-particle correlations [ 13|. It was expected that 
increasing the collision energy would increase the density and hence the expansion of the 
produced medium, so that the HBT radii should grow with y/s. 

The final interacting hadronic medium was discussed in terms of an ideal resonance 
gas, following an old suggestion [ |T^ brought to hadron physics by Hagedorn [ |^ : an 



interacting system of elementary constituents can be replaced by a non-interacting gas of 
resonances, provided the elementary interactions are resonance-dominated. This would 
provide the relative abundances of all hadron species in terms of just two parameters, the 
temperature and the baryon number density. One particularly interesting feature here 
was the fact that in elementary hadronic interactions, an overall reduction of strangeness 
production was observed. Nuclear collisions, in particular if leading to a QGP as initial 
stage were expected to remove this reduction and lead to strangenes production in 
accord with thermal predictions. 



3. Facts 

The initial energy density, as specified by the Bjorken estimate, Eq. (|I]), was measured 
in almost all SPS experiments. In Fig. ^ we show e as function of centrality, determined by 
the number of participant nucleons [ |l6|, [l^; it covers the range from somewhat above 1 to 
almost 3.5 GeV/fm^. Finite temperature lattice calculations, as already mentioned, give 
for the energy density at deconfinement, e(Tc), values around or slightly below 1 GeV/fm^ [ 
^. However, also high energy pp collisions lead to energy density estimates well above 
1 GeV/fm3. 

Thermal photon or dilepton spectra have so far not been unambiguously identified. 
Some photon excess over the expected hadronic decay yield has been reported [ |1^ , but 
its interpretation is still open. An observed excess of dileptons in the mass range between 
the and the ^/ip [ |19| has been attributed to thermal emission during the evolution of 
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Figure 2. The energy density in Pb-Pb collisions at ^/s = 17 GeV [ 16, 17 
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Figure 3. The ratio of J ftp to Drell-Yan production, measured by NA50, in Pb-Pb col- 
lisions at = 17 GeV, vs. Et (left) and Ezdc (right). The solid line indicates the 



extrapolation of the normal nuclear absorption infered from p-A collisions [ 16 



the system [ but only a small fraction would be due to the hot early stage; also here, 
more data is needed to understand the origin of the observed effects. 

J/ip production was found to be suppressed in 0-U, S-U and then Pb-Pb collisions; 
the suppression always increases with centrality [ Studying p-A collisions, it was 
observed that already normal nuclear matter leads to reduced charmonium production. 
Extrapolating this 'normal' suppression to AB interactions is enough to account for the 
observed yields up to central S-U collisions. In central Pb-Pb collisions, however, an addi- 
tional 'anomalous' suppression was observed [ 0. Peripheral Pb-Pb collisions follow the 
normal pattern; then, with increasing centrality (measured either through the associated 
transverse energy or through the number of participating nucleons), there is a pronounced 
onset of further suppression (Fig. ^). The very peripheral data points shown in this figure 
are somewhat above the normal suppression curve; this was shown to be due to beam-air 
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Figure 4. The dilepton spectrum in Pb-Au 
collisions at ^Js = 17 GeV, compared to 
the expected yield (solid line) from known 
hadronic sources [ ^ . 
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Figure 5. Hadronization temperature, T, 
and transverse flow measure, Pt, as func- 
tion of the beam energy [ . 



interactions and is no longer the case in the data collected in year 2000, when the target 
was placed in vacuum [ ^ . 

The dilepton mass spectrum in the region below the p peak was indeed found to differ 
considerably from the yield and form expected from known hadronic sources [ ^ , indicat- 



ing the presence of in-medium resonance modifications (Fig. H). This 'low mass dilepton 
enhancement' is observed in S-U and Pb-Pb collisions, and for the latter at beam energies 
of 40 GeV as well as of 160 GeV. The broadening of transverse momentum spectra, ex- 
pected as consequence of transverse flow, was observed in the predicted form of increasing 
broadening with hadron mass. If parametrized in terms of radial flow [ the extracted 
flow velocity Pt appears to saturate with increasing collision energy (Fig. ^). 

The transverse momentum spectra, moreover, also showed the azimuthal anisotropy 
predicted for non-central collisions. The behavior shown in Fig. ^ indicates at low col- 
lision energy a reduced emission in the transverse direction, where spectator nucleons 
are present; at high energy, there is enhanced production in the direction of the higher 
pressure gradient as determined by the anisotropic interaction volume [ 

The expected increase of source size in HBT correlation studies has not been observed; 
instead, one finds that, with increasing collision energy, the observed source radii are 
essentially determined by those of the involved nuclei Thus one finds -Rside — -Rout — 
5 — 6 f m for Au-Au/Pb-Pb collisions from AGS to SPS (and on to RHIC), as seen in 
Fig. 1^. This is to be contrasted to predicted values for Pb-Pb collisions ranging from 
14 fm for the SPS to 20 fm for RHIC, if freeze-out occurs when the mean free path of a 
hadron exceeds the system size [ p7 | 
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Figure 6. Elliptic flow measure, ^2, as func- Figure 7. HBT radii at different collision 
tion of the beam energy [ 25]. 



energies [ ^ . 
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Figure 8. Strange hyperon production per participant nucleon, normalized to the ratio 
from p-Be, as function of the number of participant nucleons, as measured by NA57 



For a freeze-out at the ideal pion gas energy density, one has Rf ~ {dNh/dyY^^, which 
increases from 9 to 12 fm between SPS and RHIC energies. Another unexpected feature 
is that the outward and sideward radii, i?out and i?side, whose difference could provide 
information about the life-time of the emitting medium, appear to be approximately 
equal. The HBT source radii in nuclear collisions thus seem to be energy-independent 
and fixed by the initial nuclear size [ |28[ . 

Hadron abundances showed the expected enhancement of strangeness production; in 
Fig. ^ we show the most striking example, where the production of strange baryons is 
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increased up to 10 times and more in comparison to pp rates [ In Fig. ^ it is seen 
that, quite generally, the Wroblewski measure for strangeness content is about a factor 
two higher in AA collisions than in elementary particle reactions. Of course, the role of the 
changing baryon density has to be kept in mind when comparing strangeness production 
in pp and AA collisions, or among AA collisions of different energies. 
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Figure 9. Wroblewski strangeness measure 
for different collision configurations [ ^ . 



Figure 10. Thermal hadronization and 



species abundances, from Ref. [ 31 . 
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Figure 11. Hadronization temperature as function of c.m.s. energy, for e~^e (squares), pp 
(open triangles), pp (circles) and AA collisions (filled triangles) [ |32|. 



Nevertheless, in all cases the abundance of the different hadron species is very well 



described by an ideal resonance gas. A fit to the latest NA49 data is shown in Fig. 10 [ 31 



The hadronization temperature increases with increasing collision energy and converges 
for all hadron production, from e'^e' and pp/pp to high energy A A collisions, to a stable 
value around 150-180 MeV; see Figs. ^ and [TT]. The associated baryo-chemical potential 
decreases when going from AGS to SPS and on to RHIC. 
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4. Lessons 

What can we learn from the results just listed? On the one hand, this is still a rather 
subjective issue, and so I proceed with due apologies to all those who feel that I have 
misunderstood what nature is trying to tell us. On the other hand, our theoretical under- 
standing is presently at what might be a crucial point. It was and is generally thought that 
high energy nuclear collisions will produce a thermalized medium of deconfined quarks 
and gluons, so that we can test the predictions of finite temperature statistical QCD. 
In recent years, however, the role of the primary parton state has been emphasized in 
two independent but closely related approaches, color glass condensation [ ^ and 
parton percolation [ ^ . We are thus today confronted with the interesting question 
of knowing what part of the observable features of nuclear collisions are determined by 
the primary initial conditions, before any thermalization may take place. 

Assuming thermal equilibration, with Tq — 1 fm, the initial energy densities in central 
Pb-Pb collisions at the SPS are estimated to reach 3 GeV/fm^ or more. What kind of 
medium can we have at this stage — could it be hadronic matter, a system of hadronic 
comovers of some kind? I believe that the answer is a clear NO. The average energy 
per hadron in these reactions is around 0.5 GeV, so that the mentioned energy density 
estimates would imply a hadron density of some 6 hadrons/fm^. The hadronic mean free 
path in such a medium, A = l/nhCTh, reaches something like A ~ 0.06 fm, with a typical 
hadron-hadron interaction cross section of around 30 mb. The time needed by a hadron 
to react to a single collision is around 1 fm, and in that time it travels a distance ~ 1 fm. 
In other words, the hadron undergoes almost twenty collisions in the time it needs to react 
to a single one. This rules out any medium consisting of independent hadrons [[l|, |37[| . 

If the produced medium is not hadronic, does that mean that it is a quark-gluon plasma? 
The answer is not as easy as we once thought, so let us approach this problem slowly and 
carefully. 

In a comparative study of the hadronization temperature Th as function of the initial 
energy density for different production processes, from e~^e~ and pp/pp reactions up to the 
heaviest AA collisions, we have to take into account the different values of the associated 
baryo-chemical potential. Since in all cases hadronization can be described in terms of an 
ideal resonance gas, we can extrapolate the baryon-rich low energy AA results at constant 



entropy to /i^ = 0; the resulting behavior is shown in Fig. |1^ [ |38[. It is quite striking: 
above some threshold value ec — 1 GeV/fm^, putting more energy into the system has no 
effect on the critical temperature Th which determines the hadron abundances. In the 
end, Hagedorn was right [ |^. This is like observing that water vapor always condenses 
to water at 100 °C, whatever was the initial temperature of the vapor. While this shows 
that the water at the condensation point gives no information about the earlier vapor 
temperatures, it also indicates that the condensation process exhibits critical behavior. 
In the same way, a constant Tfj for hadronic matter of arbitrarily high energy densities 
appears to signal critical behavior [ With this in mind, we now look once more at 
the experimental numbers. The saturation temperature is 150-180 MeV, and is reached 
for e ~ 1 GeV/fm^; both values are in excellent agreement with the mentioned results 
from lattice QCD studies [ For very small or vanishing baryon density, experiment 
and theory thus fully agree on the boundary of hadron physics. 



0.2 



250- 



200- 



S.150- 
E 



■B 100- 



1 50- 
X 



• e*e \'s = 91 GeV 

pp Vs = 900, 1800 GeV 
■ AA \'s = 2.6, 3.3, 4.9, 8.8, 19, 17, 130 GeV 



ni-r.i 



10 



10 



Energy density (GeV/fm ) 

Figure 12. Hadronization temperature vs. 



0.1 



T [GeV] 


• pp/pp 
T SPS 
■ AGS 
\ ♦ SIS 


thermal \ \ 
freeze-out\ \ 

\l 


\ dcconfincincnt 

chemical \ 
freeze-out \ 


nuclear matter 

\ 


\ |i [GeV] 

L) 1 \ 1 


0.5 


1.0 1.5 



the initial energy density, from Ref. [ 38 . 



Figure 13. Deconfinement, chemical and 
thermal freeze-out vs. T and /i [ HTI, . 



Extending the resonance gas analysis of species abundances to finite baryon density, 
one obtains the 'chemical' freeze-out curve [ ^] in the T — fis plane; it is illustrated 
schematically in Fig. |T3|. Also included in this figure are the thermal freeze-out curve 
(obtained from a flow analysis of transverse momentum spectra) and the expected decon- 
finement transition curve [ ^ . Some first and quite promising attempts of lattice gauge 
theory to study deconfinement at non-zero baryo-chemical potential, fiB, agree with the 
deconfinement curve for fiB < 0.8 GeV [ It is clear that the deconfinement curve 

must lie above the freeze-out curve at large /i^: the interactions in cold dense nuclear mat- 
ter are dominated by baryon repulsion and hence cannot be described in terms of an ideal 
resonance gas. At T = 0, freeze-out occurs at normal nuclear matter density, where baryon 
repulsion and attractive (pion exchange) interactions just balance each other. Neutron 
stars are examples of confined states of higher baryon densities. 



Since Figs. 0, |T2| and [T^ include both data from elementary particle and from nuclear 
collisions, statistical hadron abundances alone are not sufficient to indicate the presence 
of a large-scale system. Therefore, it is necessary to check if AA collisions produce large 
thermal systems or just the sum of the smaller ones. The answer to that question is 
obtained through the study of strangeness production. In elementary particle collisions, 
there is an overall suppression of strange particle production compared to the predictions 
of the grand canonical form of an ideal resonance gas. This can be interpreted as the 
consequence of local strangeness conservation: if only one ss pair is formed in the gas, the 
corresponding Boltzmann factor should be exp{— 2ms/T}, instead of the grand canonical 
form exp{— ms/T}. If nuclear collisions lead to a connected large-scale system, strangeness 
conservation can occur between the products from different nuclei, leading to the grand 
canonical Boltzmann factor. As a consequence, we expect for the Wroblewski measure in 
elementary particle collisions 

2 < ss > 

= ^ -^^^WJ^ ^ exp{-2m,/T} ^ 0.2, (3) 
< uu > + < dd> 

while central high energy AA collisions should approach 

A, ~ exp{-m,/T} ~ 0.4. (4) 



Keeping in mind that the different baryon densities also lead to modifications, we see in 
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Fig. I that Eqs. (|]) and (|) are in fact quite well satisfied. We thus conclude that nuclear 
collisions indeed produce large-scale connected systems. 

The observed behavior of the hadronic transverse momentum spectra is also often taken 
as indication for collective behavior of the produced medium. In particular, azimuthal 
asymmetries seem to indicate that global effects influence the spatial hadron emission: in 
non-central collisions, there is more (elliptic) flow in the direction of the greater pressure 
gradient. For the time being, however, it is not so easy to reconcile a constant radial flow 
velocity at high energy (Fig. |^) and energy-independent HBT radii of nuclear size (Fig. |^) 
with an initial state of partonic matter at higher and higher energy densities. Apparently 
a change of conditions in the pre-hadronic phase, whatever its nature, is not reflected in 
radial flow or source size. 

We conclude that there seems to be a universal limit to the hadronization regime in 
high energy nuclear collisions; the abundances of the produced hadron species, including 
the strange ones, are determined by values of temperature and energy density in good 
agreement with what is known about the confined region of the QCD phase diagram. 
As function of the incident energy, the hadronization temperature Th, the strangeness 
content and the transverse momentum broadening (radial flow velocity (3t) all seem to 
converge to constant values. The value of the temperature Th — 150 — 180 MeV and of 
the corresponding energy densities agree with the confinement /deconfinement transition 
values obtained in lattice QCD. The strangeness content agrees with what is expected 
from the value of the Boltzmann factor exp{— m^/T} at T = Th- In fact, even the 
low mass dilepton enhancement appears to remain essentially invariant under changes of 
nuclear size and of collision energy (in the SPS range). The origin of the effect is not yet 
really clear (see, however, [0), but low mass dilepton production does seem to behave in 
the same way as the mentioned hadronic observables. Since it may be the most sensitive 
to changes in the initial baryon density [|4^, this observation could be a reflection of the 



weak /i-dependence of the confinement limit in the small fx region. 

Before analyzing the earlier stages of the nuclear collisions, we briefly mention attempts 
to include charmonium production in a statistical hadronization scheme, either by assum- 
ing thermal formation at the freeze-out point T = Th of the resonance gas [ ^ , or by 
assuming a thermal distribution of the produced cc pairs among the different possible 



open and hidden charm channels [ To evaluate the viability of these approaches, we 
recall that J/ip production from pp to central S-U and peripheral Pb-Pb collisions is seen 
to be in accord with the normal nuclear suppression determined in p-A interactions. In 
contrast, more central Pb-Pb collisions show a well-defined onset of anomalous suppres- 
sion in a narrow centrality range, seen both in the ratio of J/ip to Drell-Yan production 
and in the overall J/ip yield (see Fig. Such an onset cannot be obtained in any thermal 
production model, nor can the centrality-dependent ip'/ip ratio found in S-U collisions [ 
Hence, the present data indicate that the production and subsequent fate of charmo- 
nium states in nuclear collisions differ significantly from that of light hadrons. We shall 
return shortly to the interpretation of the observed charmonium patterns. 

The study of the behavior of (light quark) hadrons has thus provided us with clear 
indications determining where the hadronic regime ends. It does not tell us, however, 
where the QGP begins. High energy elementary particle collisions lead to initial energy 
density estimates higher than those of Pb-Pb collisions at SPS or RHIC energies (see 
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Figure 14. Four-jet and two-jet decays of W^W pairs in e^e annihilation at LEP. 



Fig. |T2D; for the resulting hadronization temperatures we always get the universal T^-. 
Therefore, one has tried to find evidence for pre-hadronic interactions, between partons 
from different sources. In e~^e~ interactions at LEP energies {^/s ~ 200 GeV), this can be 
done by comparing hadron multiplicities in two-jet and four-jet decays of W^W~ pairs 
produced essentially at rest (see Fig. |n|). 'Cross talk' between jets from different W^s 
should be reflected in hadron multiplicities, with the prediction [ Q that Nh{Ajet) < 
2Nh{2jet). The data do not show any evidence for such 'color interconnection' even 
though the energy density calculated according to the Bjorken estimate (|lD reaches values 
around 2 GeV/fm'^. The partons from different sources simply don't seem to know of each 
other. Color interconnection is clearly a prerequisite for thermalization on a partonic level; 
hence e > 1 GeV/fm^ is a necessary, but not a sufficient condition for QGP formation. We 
must somehow assure that the partons in the early stages of the collision reach a density 
so high that they interact and allow color interconnection to set in. In other words, the 
pertinent question is: when do the individual partons from the different nucleons of the 
colliding nuclei condense to form a large connected cluster? 



From percolation theory [ it is known that the formation of large-scale clusters 
is a critical phenomenon; it does not occur gradually as function of parton density. In 
the 'thermodynamic' limit of large nuclei, the cluster size diverges at a specific critical 
parton density, and even for finite systems it suddenly increases in a very narrow band 
of density. Thus there exists a statistical approach to critical behavior which does not 
pre-suppose thermalization and which can be applied to the pre-equilibrium stage of the 



nuclear collision evolution \ 35, 36 



In a central high energy nucleus-nucleus collision, the partons in the two nuclei lead in 



the transverse plane to an initial condition schematically illustrated in Fig. 15. The trans- 





Figure 15. Parton distribution in the transverse plane of a nucleus-nucleus collision, for 
small (left) and large (right) nuclei. 
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verse size of the partons is essentially determined by the intrinsic transverse momentum, 
cTq ~ Tt/k'^. The number of partons contained in a nucleon is known from deep inelastic 
scattering experiments. It is parametrized by a parton distribution function depending 
on the fraction x = k/p (denoting parton and nucleon momenta by k and p, respectively) 
and on the scale Q used to resolve the nucleonic parton structure. While in lepton-hadron 
scattering the scale is set by the virtual photon, in minimum bias nucleon-nucleon or AA 
collisions it is determined by the transverse momentum of the partons themselves. 

We denote the nuclear radius in Fig. [1^ by R, the average parton radius by r, and 
then study the variation of the average cluster size as function of the parton density 



n 



N, 



parton 



/ttR . In the limit i? — > cxd, the cluster size diverges at the percolation 



threshold Up ^ l.lS/vrr 



Sr. 



n] 



(5) 



with the critical exponent 7 = 43/18 [ Percolation thus specifies the onset of con- 



nection as a critical phenomenon. For finite i?, there is a pronounced but finite peak at 



a slightly shifted density, as shown in Fig. 16 
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Figure 16. Average cluster size Sd as function of parton density n; np denotes the perco- 
lation point in the limit R ^ 00. 



To apply this formalism to nuclear collisions [ 36 1, we need the transverse parton size 
and the effective number of partons for a given collision configuration. The distribution 
of partons in a nucleon is determined in the analysis of deep inelastic scattering data; per 
nucleon in a nucleon-nucleon collision, at mid-rapidity and for c.m.s. energy ^/s, one has 



parton 



. dy 



(6) 



y=0 



where g, q, q label the gluon, quark and antiquark distributions, respectively, and the 
sum runs over the quark species. At y = 0, the fractional parton momentum x = k/p 
becomes x = Qj \fs. As shown in Fig. |15|, there is a distribution of partons of different 
transverse sizes up to the resolution scale Q^; we approximate this by the average value 
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In nucleus-nucleus collisions at SPS energies, it is a good approximation to consider 
that the activated partons originate from wounded or participant nucleons. At higher 
energies, there will also exist collision-dependent contributions [ Q . 

We thus obtain for central AA collisions the percolation condition 



it determines the onset of color connection as function of A and of the cm. collision energy 



The parton condensate formed beyond the percolation point consists of overlapping and 
hence interacting partons from all involved nucleons. It thus constitutes a deconfined pre- 
thermal medium, which is a necessary precursor for QGP formation, since color connection 
is a prerequisite for thermalization. The parton condensate is characterized by a scale 
Q = Qs, where Qs{A, i/i) specifies the onset of percolation for central collisions at a given 
A and y/s. The average transverse momentum of partons in the condensate is determined 
by this scale, so that Qs is in a sense a precursor of temperature. Increasing A or ^/s 
leads to higher Qg, indicating something like a "hotter" medium. 

We have here concentrated on the onset of parton condensation as a critical phenomenon 
determined by percolation of quarks and gluons as geometric entities in the transverse 
plane. The behavior of the dense parton condensate in the limit of large A and a/s, the 
color glass condensate, has been studied intensively over the past years in terms of classical 
color fields. It is becoming increasingly clear from these studies that the pre-equilibrium 
stage plays a much more decisive role in nuclear collisions than previously envisioned. In 
the last part of this section, we consider one interesting experimental consequence of the 
onset of parton condensation. 

To study this onset as function of centrality in a given A — A collision, we have to 
replace 2A/'7iR\ in Eq. (|^ by the density of wounded nucleons n^^b) calculated for a 
given impact parameter b, using the actual nuclear profile (Woods-Saxon). For Pb-Pb 
collisions at i/i =17 GeV, this leads to 



Instead of the impact parameter, we have here used the number of wounded nucleons to 
specify the collision centrality, with N^, — 390 for 6 = 0. Note that at this threshold point, 
the Bjorken estimate for the energy density gives ~ 2.5 GeV/fm^. This clearly shows that 
the onset of color connection only occurs at a much higher energy density than the value 
given by finite temperature lattice QCD for a thermalized medium. 

As mentioned above, charmonium suppression in nuclear coUisons was one of the pro- 
posed signatures for quark-gluon plasma formation. In the light of present thinking, we 
have to consider the fate of the charmonium states already in the parton condensate 
phase. This can be addressed in different ways, invoking generalized color screening [ 
dipole break-up in a random color field [ etc. We shall here simply consider the 
intrinsic scale Qi of a given charmonium state i (J/V', Xc, 4'') assume that when it is 




(7) 



v/i. At = 20 GeV, we obtain [ 



for y4 < 60, no color connection, 

for A > 60, a connected parton condensate. 



for A^^ < 150, no color connection, 

for A^^ > 150, a connected parton condensate. 
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exceeded by the condensate scale, 

Qs > Q^, (8) 

the charmonium state is dissolved. We are thus assuming that plays in the pre- 
equilibrium state the role of the critical temperature or of the screening mass in a thermal 
medium. 

To show the consequences of this assumption, we must first recall that J/ip production 
in hadronic collisions occurs in part through feed-down from higher excited states: about 
60% of the observed J/ip in proton-proton collisions are directly produced (IS) states, 
the remaining 40 % coming from Xc (~ 30 %) and tp' (~ 10 %) decays. The intrinsic scales 
of these states are given by their radii. 



rj/v, ~ (0.9 GeV)- 



'x — 



(0.6 GeV)- 



(0.6 GeY)-\ 



(9) 




100 150 200 250 300 
Number of participants 



400 



Figure 17. Centrality dependence of the percolation scale, Qs, in Pb-Pb collisions, at SPS 
and RHIC energies [ p6| . 



From the centrality dependence of Qs shown in Fig. |T^, we find that, at SPS energies, 
charmonia suppression sets in at 

^ 150 for the Xc and ip' contributions, 
Nw ^ 250 for direct J/ip production, 

to be compared with the observed pattern shown in Fig. |18[ Obviously these results must 
be studied in more detail, but it appears quite likely that the onset of J/ip suppression 
indeed indicates an onset of deconfinement, while not implying any thermalization. 



5. Summary 

So let us summarize what we have learned in the first fifteen years of ultra-relativistic 
heavy ion collisions. 

We know that there is an intrinsic limit to hadronic matter, in accord with the phase 
diagram determined in statistical QCD. The confinement boundary is established by the 
SPS/AGS program and agrees, both quahtatively and quantitatively, with the predictions 
of lattice QCD. I believe that the RHIC and LHC experiments can only reconfirm this. For 
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Figure 18. ^/if) production in pp, p-A and AB collisions, as function of the number of 
wounded nucleons [ [If 



small we know the onset of hadronization, and we know that strangeness enhancement, 
transverse momentum broadening and HBT radii saturate at this onset. 

The early stage of the medium produced in nuclear collisions at the SPS is partonic. 
With increasing A, the partons begin to form connected clusters, and at a certain critical 
density, they form a color condensate consisting of interacting partons from many different 
nucleons. Hence, at this point, color deconfinement begins: partons no longer have a clear 
origin or exist in well-defined numbers. The parton condensate is a necessary precursor 
of the QGP: it provides the interacting deconfined partons which, if given enough time 
for equilibration, can make a QGP. 

Charmonium states of different binding energies and spatial sizes probe different par- 
ton scales of the produced initial state. The observed step-wise form of anomalous ^/ip 
suppression appears to provide the first signal for color deconfinement through parton 
condensation; this does not require any thermalization. The forthcoming ^/ip production 
measurements at full SPS energy but with lower A, to be performed by NA60, should 
further clarify this point. In particular, the comparison between Pb-Pb and In-In ^/ip 
suppression patterns should determine the critical densities and scales. 

Pioneering results in physics always need to be reconfirmed. But if these limits of 
confinement survive the test of time, we will remember that they were first found in 
SPS /AGS experiments. 

So, at the beginning of the collision evolution at the SPS, there appears to be color 
deconfinement; at the end, thermalization and collective behavior. To produce a quark- 
gluon plasma medium, we need to have both at the early stage. Let us see if the future 
data from experiments at the much higher RHIC and LHC energies can achieve this. But 
no matter how well-defined the road for further exploration may seem to be, I am sure 
that the forthcoming studies will rediscover the one feature which has made this field so 
challenging and exciting. It is perhaps best summarized by the Spanish poet Antonio 
Machado [^■. 

Caminante, son tus huellas el camino, y nada mas; 
caminante, no hay camino, se hace camino al andar. 

Traveller, the road is nothing more than your footprints; 
traveller, there is no road, you make it as you go. 
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